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Abstract 

A new method is described for determining the quantum correlations at differ- 
ent times in optical pulses by using balanced homodyne detection. The signal 
pulse and sequences of ultrashort test pulses are superimposed, where for chosen 
distances between the test pulses their relative phases and intensities are varied 
from measurement to measurement. The correlation statistics of the signal pulse 
is obtained from the time-integrated difference photocurrents measured. 
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Methods for measuring the quantum statistics of optical fields have been of 
increasing interest. One of the most fruitful experimental methods has been 
balanced homodyne detection JIJ including homodyne correlation measurements 
[0. In particular, combining a single-mode signal field and a strong local oscillator 
(LO) by a 50%:50% beam splitter and measuring the interfering fields in the two 
output channels by means of two photodetectors, the difference- count statistics 
yields directly the field-strength statistics of the signal mode for a certain phase 
parameter. Since knowledge of the probability distributions of the field strengths 
for all phases within a 7r-interval is equivalent to knowledge of the quantum state 
PI, the field-strength distributions can be used to reconstruct the quantum state 
in terms of the density matrix or other representative quantities, such as phase- 
space functions. 

In the pioneering experiments by Smithey, Beck, Raymer, and Faridani the 
Wigner function was determined from the measured field-strength distributions 
by means of inverse Radon transform and the density-matrix elements in a field- 
strength basis were then obtained as Fourier transforms of the Wigner function 
(see also Ref. 0). The method also called optical homodyne tomography (OHT) 
has been improved in order to obtain the density matrix in a more direct way 
0. In particular, the problem of direct sampling of the density matrix in the 
photon-number basis has been studied in a number of papers 0. The method 
can be extended to the detection of multimode light by increasing the number of 
the input (and output) channels of the apparatus. The problem of reconstruction 
of the quantum state of two correlated modes was first considered in Ref. ||, and 
an analysis of the determination of the quantum state of a correlated iV-mode 
field by multiport homodyning is given in Refs. ||. 

In the experiments reported in Ref. [|J optical pulses are used and the time- 
integrated difference-count statistics is measured. Hence, a single-mode density 
matrix is reconstructed from the measured data, which only yields information 
on the quantum statistics of the signal pulse as a whole, but does not contain 
any information about its internal quantum statistics. Recently, Munroe, Bog- 
gavarapu, Anderson, and Raymer have performed OHT experiments using 



LO pulses that are short compared with the signal pulses. In this way they have 
been able to obtain the photon-number statistics of a signal pulse at different 
times in the pulse. 

In this paper we propose a method for measuring also the internal correlation 
statistics of optical pulses, which applies to both classical and nonclassical light. 
In particular, the method can be used to study correlation-assisted coherence 
properties in an experiment analogous to Young's interference experiment, but 
with a "temporal double slit" inside a pulse. Another example is the measurement 
of the photon-number correlation at different times in the pulse. It is worth 
noting that even in classical optics the measurement of the internal amplitude- 
phase structure of pulses is not trivial and has been a subject of current interest 
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To measure the internal quantum statistics of an optical pulse, an apparatus 
is desired that performs an appropriate mode decomposition of the pulse with- 
out introduction of additional noise and measures the quantum statistics of the 
correlated modes. This can be achieved by combining the signal pulse, through a 
50%:50% beam splitter, with a train of N LO pulses that are short compared with 
the signal pulse and measuring the time-integrated difference-count statistics in 
the two output channels of the beam splitter (Fig.[|). The signal pulse consists, 
in general, of a continuum of monochromatic modes with photon destruction and 
creation operators b(u>) and b>(u)'), respectively, [b(u),b'(u')] = 5(u' — u). Let us 
suppose that the strong LO pulses are prepared in coherent states, so that the 
positive-frequency part of the kth pulse centered at time t k can be described by 
a function r Yk9k{t)i where gk{t) is assumed to be normalized to unity and j k is a 
complex number. We now consider the operators 

a k = J dug* k (u)b(u), (1) 

where gu{oS) is the Fourier transform of gk(t), gk(uj)= (2tt)~ 1 ^ 2 J dio gk(t) exp(iut). 
It can be shown that when the LO pulses guit) and gk'(t), k ^ k', are (ap- 
proximately) nonoverlapping, then the operators a k and a k , satisfy the standard 
bosonic commutation relation [a^a^/] = 5 kk i. Hence, the train of LO pulses can 



be used to introduce a set of nonmonochromatic modes [12] and "probe" the 
correlation statistics of the signal pulse in terms of these modes. In particu- 
lar, measurement of the integrated difference- count statistics can be shown to be 
equivalent to measurement of the sum of signal-pulse field strengths 

F = "£q k F k (tp k ). (2) 

k 

Here, 

F k {<p k ) = 2- 1 / 2 (a k e~^ + a\e^) (3) 

is the field strength associated with the kth pulse-like nonmonochromatic mode 
centered at time t k in the signal pulse and q k is a non-negative real parameter 
describing the relative (square root of the) intensity of the kth LO pulse. Mea- 
suring the distribution of the sum field strength (|3]) for a sufficiently large set 
of values of <p k and q k , we may obtain the quantum statistics of the signal pulse 
within the frame of an iV-mode density matrix. 

The train of LO pulses may be obtained from one generating LO pulse using 
interferometric methods, so that the distances between the LO pulses, their inten- 
sities and relative phases can be controlled. Let us further assume that the signal 
pulse can be appropriately triggered by the generating LO. In this way, the times 
t k in the signal pulse and the absolute values and the relative phases of the F k 
can be controlled. When the "positions" t k of the LO pulses are shifted towards 
t k + At, then the phases Lp k are shifted towards ip k + uoAt, ujq being the carrier 
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frequency of the signal and LO pulses. Since times of the order of magnitude 
of Uq 1 cannot be controlled, an absolute phase control seems to be impossible 
when the signal and LO pulses come from different sources. In this case, only 
quantities averaged with respect to the unknown phase can be obtained. 

Let us consider the simplest case when the signal pulse and two short nonover- 
lapping LO pulses centered at the times t\ and t 2 in the signal pulse are superim- 
posed. In this case, measurement of the time-integrated difference-count statistics 
yields the statistics of a sum of two-mode field strengths of the signal pulse, 

F = F(<p,Av,q) = F 1 (v) + qF 2 ( V + A<p). (4) 

The measured moments of F are related to the moments and correlation functions 
of the signal pulse at the two times t± and t 2 in the pulse as 



( Fn ) = £ [ k )l*( F i(<p)F 2 *(<P + Ap)). (5) 

Measuring (F n ) for n + 1 values of q, from Eq. (^) one obtains a set of n + 1 
linear algebraic equations whose solution yields the signal-pulse moments and 
correlation functions (F™~ k (ip)F 2 (ip + Aip)) for the chosen values of cp and Aip. 
Varying tp and A<p, the procedure can be repeated many times to obtain the (p 
and Aip dependences of the moments and correlation functions. In the limit when 
the number of measurements goes to infinity all the moments and correlations 
can be obtained, the knowledge of which is equivalent to knowledge of the two- 
mode density matrix. Clearly, when the phase ip cannot be controlled the phase- 
averaged (even) moments 

W) = ±-Jd<p (F n ) (6) 
can only be measured and Eq.© modifies to 



n being even. 

Let us briefly discuss the situation of nonperfect detection. In this case, ak 



must be replaced by rjau + Jr](l—r])ck (cf. Ref. |0|)> where r\ is the detection 
efficiency (rj < 1) and dk and c\ are bosonic noise operators. The corresponding 
noise modes can be considered to be in vacuum states, so that for each of these 
modes the mth moment of the corresponding field strength is given by 2 _m//2 (m — 
1)!! for even m and zero for odd m [{m— 1)!! = 1 • 3 • • • (m — 3)(m — 1) for m > 2; 
(— 1)!! = 1]. We then find that in the sum in Eq. ([?]) the replacement 



(Fr%)F 2 % + A<^)) 




x (n-k-l-l)\\ {k-m-l)\\ {Fl(ip)F 2 m (tp + Aip)) (8) 
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must be made. When k is even (odd) the I- and m-sums in Eq. (|j) run over 
the even (odd) integers between zero and n — k and between zero and k, respec- 
tively. The replacement means that the reconstructed moments and correlations 
in Eq. fl7|) are attenuated and "contaminated" by the moments of lower orders. 
Nevertheless, we can recurrently correct them by subtracting the contaminating 
terms, because only lower-order terms are included. 

To illustrate the method, let us briefly consider the phase-averaged second- 
and forth-order moments of the measured time-integrated difference- count statis- 
tics. Information about the lowest-order correlation in the signal pulse at two 
times t\ and t 2 is obtained from the second-order moment, 



= (FKv)) + q 2 (F$(ip)) + 2q(F 1 (i P )F 2 (<p + Aip)). (9) 



Whereas (Fi(ip)) and (F 2 (p)) are closely related to the mean numbers of photons 
at the two times in the signal pulse 



(*2(vO> = <**>+ 5 (10) 

(k= 1,2), the cross term 



{F l {i P )F 2 {i P + Aip)) = ±(a\a 2 )e- iA v + c.c. (11) 

as a function of the difference phase A<p can be regarded as the second-order co- 
herence function that probes the effect of signal-pulse interference at a "temporal 



double slit" in the pulse [14] . Choosing q = 1 and measuring the second-order 

moment (F 2 ) for a set of phase differences Aip and Aip + 7r, the interference 
term is simply given by the difference of the measured moments at the two phase 
differences. In the case of nonperfect detection the photon-number terms are 
attenuated and contaminated by a constant, 



(m^^v'iF^ + H^-v), (12) 
whereas the interference term is only attenuated, 



(F 1 (^)F 2 ( ¥ > + A<p)) -> ri 2 (F 1 {ip)F 2 { ( p + Aip)). (13) 



Extending the measurement to the fourth-order moments (F 4 ) for different 
values of q and Aip, the phase-averaged fourth-order moments and and correla- 
tions of F 1 ( ( / 3 ) an d F 2 (ip>) can be obtained. The moments (Ff{tp)) and (F 2 (p)) 
carry information about the photon-number statistics (first- and second-order 
moments of the photon numbers) of the signal pulse at two different times in the 
pulse. The two-time correlations are given by 



(F?(<p)F 2 (<p + Aip)) = I (F^F^tp + Aip)) + |((a t 1 2 a 1 a 2 )e-^ + c.c), (14) 



4 



(F 1 ^)Fi(i P + A<p)) = | {F 1 {i P )F 2 (i P + A<p)) + |((a t 2 2 a 2 a 1 )e iA ^ + c.c.) (15) 



with (F 1 ((p)F 2 ((p + Aip)) from Eq. (O), and 



(F?(ip)Fl(ip + A(p)) = \ ((a\ 2 al)e-^ + c.c.) + (fa + |)(n 2 + ±)>. (16) 



In particular, from Eq. (16) together with the mean photon numbers obtained 
from the second-order moments, Eq. (|Il]), the photon- number correlation at the 
times t\ and t 2 in the signal pulse can be determined, which for stationary ergodic 
fields is usually measured in the Hanbury-Brown and Twiss experiment. When 
the detection efficiency is less than unity the terms (F^(if)) changes as 



<^%)> - V'(F^)) + 3r/ 3 (l - V)(F^)) + h 2 (l - V) 2 , (17) 



and with respect to the symmetrical cross term (|I6) we find that 



(FHV)F 2 2 (V + Av>)> -> V A (F 2 ((p )F 2 2 (ip + Ay?)) 

+^ 3 (1 - rj)[(t?(<p)) + (F$(cp))} + ^ 2 (1 - (18) 

The inclusion in the analysis of higher-order moments is straightforward. In 
summary, the proposed method may serve as a useful tool for determining the 
correlation statistics of optical pulses, with special emphasis on nonclassical light. 
The application to pulses becomes feasible since only ensemble averaging is used. 
The other advantage of the method is that it also applies to weak quantum light, 
because homodyning (with strong LO) "amplifies" the signal. In particular, the 
method can be used to obtain the correlation statistics of pulses of the type used 



in Ref. JT0[ for studying the photon-number statistics at different times in a pulse. 
In the paper we have restricted attention to the determination of moments and 
correlations at different times in the signal pulse. The problem of reconstruction 
of the associated density matrices will be considered in a forthcoming paper. 
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Figure 1: Scheme of the measurement. The signal pulse and a train of strong 
local-oscillator (LO) pulses that are short compared with the signal pulse are 
combined by a 50%:50% beam splitter, and the time-integrated difference-count 
statistics in the output channels is measured. The train of LO pulses is produced 
interferometrically, so that the pulse distances, relative phases, and intensities 
can be controlled. 
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